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ABSTRACT 

We present results from a very deep (650 ks) Chandra X-ray observation of the galaxy group 
NGC 5813, the deepest Chandra observation of a galaxy group to date. Earlier observations showed 
two pairs of cavities distributed roughly collinearly, with each pair associated with an elliptical shock 
front. The new observations confirm a third pair of outer cavities, collinear with the other pairs, 
and reveal an associated outer outburst shock, with measured temperature jump, at ~ 30 kpc. This 
system is therefore unique in exhibiting three cavity pairs, each associated with an unambiguous shock 
front, arising from three distinct outbursts of the central AGN. As such, it is particularly well-suited 
to the study of ongoing AGN feedback. The implied mean kinetic power is roughly the same for 
each outburst, demonstrating that the average AGN kinetic luminosity can remain stable over long 
timescales 50 Myr). The two older outbursts have larger, roughly equal total energies as com¬ 
pared with the youngest outburst, implying that the youngest outburst is ongoing. We find that the 
radiative cooling rate and the mean shock heating rate of the gas are well balanced at each shock 
front, suggesting that AGN outburst shock heating alone is sufficient to offset cooling and establish 
AGN/ICM feedback within at least the central 30 kpc. This heating takes place roughly isotropically 
and most strongly at small radii, as is required for feedback to operate. We suggest that shock heating 
may play a significant role in AGN feedback at smaller radii in other systems, where weak shocks are 
more difficult to detect. We find non-zero shock front widths that are too large to be explained by 
particle diffusion. Instead, all measured widths are consistent with shock broadening due to propaga¬ 
tion through a turbulent ICM with a mean turbulent speed of ^70 km s -1 . Finally, we place lower 
limits on the temperature of any volume-filling thermal gas within the cavities that would balance the 
internal cavity pressure with the external ICM. The most stringent limit we find is kT > 16 keV. 

Subject headings: galaxies: active — galaxies: clusters: general — galaxies: groups: individual 
(NGC5813) — galaxies: individual (NGC5813) — X-rays: galaxies 


1. INTRODUCTION 

Early Chandra and XMM-Newton X-ray observations 
revealed that the amount of gas cooling to very low tem¬ 
peratures at the centers of cool core clusters is much 
less than what is expected from simple radiative cooling 
models (David et al. 2001; Peterson et al. 2001; Peterson 
& Fabian 2006). The implication is that the diffuse X- 
ray emitting gas must be heated, either by pre-heating 
during cluster formation or by ongoing energy injection. 
The most likely heating mechanism is feedback due to 
energy injection by the central active galactic nucleus 
(AGN) of the cD galaxy (see McNamara & Nulsen 2007 
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and Fabian 2012 for recent reviews). During this pro¬ 
cess, matter is accreted onto the central supermassive 
black hole (SMBH), which drives powerful jets. These 
jets evacuate cavities in the intracluster medium (ICM), 
which can drive shocks as they are inflated and subse¬ 
quently rise buoyantly (Churazov et al. 2001). The en¬ 
ergy contained in cavities and shocks is then available to 
heat the ICM, which lowers the cooling rate of the gas 
and subsequently the SMBH accretion rate. The ensu¬ 
ing decrease in AGN heating allows the gas to once again 
cool and accrete onto the SMBH, establishing a feedback 
loop that regulates the temperature of the ICM. Several 
studies have shown that, generally, the total enthalpy in 
cavities in cool core systems is sufficient to offset radia¬ 
tive cooling in individual galaxies, galaxy groups, and 
clusters (Birzan et al. 2004; Rafferty et al. 2006; Nulsen 
et al. 2007; Hlavacek-Larrondo et al. 2012). However, 
the details of how and where this energy gets transferred 
to the ICM are unclear. Weak AGN outburst shocks are 
also expected to heat the ICM, although they are difficult 
to detect and clear examples are very rare. 

Galaxy groups provide an excellent opportunity to 
study AGN feedback. Their lower temperatures are more 
easily measured with modern high angular resolution X- 
ray satellites, and their central AGN can more easily dis¬ 
turb the diffuse gas due to their shallower gravitational 
potentials as compared with clusters. Here we report on 
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results from a very deep Chandra observation of the cen¬ 
tral galaxy in the galaxy group NGC 5813 (N5813). This 
group is a relatively isolated sub-group in the NGC 5846 
group (Mahdavi et al. 2005; Machacek et al. 2011). In 
Randall et al. (2011, hereafter Rll) we presented results 
based on an initial 150 ks Chandra observation of N5813. 
The ICM in this group has a remarkably regular mor¬ 
phology, with three pairs of roughly collinear cavities and 
associated surface brightness edges, and shows no clear 
signs of a recent merger event. (Note that we refer to 
the diffuse gas in this group as the ICM rather than the 
intragroup medium (IGM) throughout to stress the con¬ 
nection with feedback in clusters and to avoid confusion 
with the intergalactic medium.) With clear, cleanly sepa¬ 
rated signatures from three distinct outbursts of the cen¬ 
tral AGN and no other significant dynamical processes 
at work, N5813 is uniquely well-suited to the study of 
AGN feedback. In this work, we focus on the implica¬ 
tions for AGN feedback and the outburst history of the 
central SMBH. 

We assume an angular diameter distance to N5813 of 
32.2 Mpc (Tonry et al. 2001), which gives a scale of 
0.15 kpc/". All uncertainty ranges are 68% confidence 
intervals (i.e., lcr), unless otherwise stated. 

2. CHANDRA DATA ANALYSIS 
2.1. Observations and Data Reduction 

The Chandra observations that were used in the analy¬ 
sis we present here are summarized in Table]!] The aim- 
point was on the back-side illuminated ACIS-S3 CCD 
for each observation. All data were reprocessed from 
the level 1 event files using Cl AO and CALDB 4.4.3. 
CTI and time-dependent gain corrections were applied. 
LC-CLEAN was used to remove background flare^J The 
mean event rate was calculated from a source free region 
using time bins within 3a of the overall mean, and bins 
outside a factor of 1.2 of this mean were discarded. There 
were no periods of strong background flares. The cleaned 
exposure times are given in Table [l] for a total time of 
635 ks. 

Diffuse emission from N5813 fills the image FOV for 
each observation. We therefore used the caldbE3 
blank sky background files appropriate for each obser¬ 
vation, normalized to match the 10-12 keV count rate in 
our observations to account for variations in the parti¬ 
cle background. To generate exposure maps, we used a 
MEKAL model with kT = 0.7 keV, Galactic absorption, 
and abundance of 30% solar at a redshift z = 0.006578. 

3. IMAGE ANALYSIS 

The exposure corrected, background subtracted, 0.3- 
3 keV Chandra image of the central region is shown in 
Figure [l] To enhance the visibility of the diffuse emis¬ 
sion bright point sources were removed, and the regions 
containing point sources were “filled in” using a Poisson 
distribution whose mean was equal to that of a local an¬ 
nular background region. To better show faint structure, 
particularly in the outer regions, we fitted the X-ray im¬ 
age with a 2D beta-model using the software package 
SHERPA (Freeman et al. 2001) and divided the image by 
the model to produce a residual image. The smoothed, 

9 http://asc.harvard.edu/contrib/maxim/acisbg/ 

10 http://cxc.harvard.edu/caldb/ 


point source free image and the residual image are shown 
in Figure [2] 

These images clearly show three pairs of roughly 
collinear cavities, along an axis from the SW to the NE. 
The inner cavity pair is surrounded by bright rims (at 
about 1 kpc), and each of the intermediate and outer cav¬ 
ity pairs is associated with an elliptical surface brightness 
edge (at ~10 kpc and ~30 kpc, respectively). In Rll, 
we showed that the 1 kpc and 10 kpc edges are shock 
fronts that were driven during the expansion phase of 
their associated cavities. The new, deeper observations 
clearly show a third (outer) cavity pair and edge, only 
hinted at in earlier observations, and allow us to unam- 
biguousl y iden tify the outer edge as an associated shock 
front (§ |4.2.2|) . These outer features are most clearly 
seen in the residual image (Figure [2]. Thus, these deep 
X-ray observations show clear signatures from three dis¬ 
tinct outbursts of the central AGN. We note that, here 
and throughout, by “outburst” we refer to the creation 
of a cavity pair and its associated shocks, rather than a 
rapid increase in AGN jet power, since these features are 
in principle con siste nt with either a constant or variable 
jet power (see § 6.2). 

To look for structure in the faint emission at large 
radii, beyond the FOV of Chandra , we examined the 
available archival XMM-Newton observations of N5813. 
N5813 was observed with XMM-Newton on three occa¬ 
sions: July 23 2005 (ObsID 0302460101), February 11 
2009 (ObsID 0554680201), and February 17 2009 (Ob¬ 
sID 0554680301). The event files were calibrated using 
the XMM-Newton Science Analysis System (SAS) ver¬ 
sion 13.5.0, and the most recent calibration files as of 
2014 July. The calibrated, cleaned event files were pro¬ 
duced after filtering for particle background flares. This 
resulted in cleaned exposure times of 120 ks, 127 ks, and 
83 ks for the MOS1, MOS2, and PN detectors, respec¬ 
tively. Details of the data reduction are described in 
Bulbul et al. (2012a). 

The background subtracted, exposure corrected, 
merged XMM-Newton image is shown in Figure [ 3 ] beside 
the Chandra image. The 10 kpc shock edges are visible, 
as are the intermediate cavities and the NE outer cavity. 
The inner ^ 1 kpc cavities and rims are not resolved. A 
smoothed version of this image, with the intensity scale 
chosen to better show faint emission at large radii, is 
shown in Figure [4j The extended emission is roughly az- 
imuthally symmetric, and shows no clear evidence of a 
fourth pair of cavities or shock edge beyond the 30 kpc 
edge detected with Chandra. Since, for the purposes of 
this study, we are most interested in the detailed struc¬ 
ture in the ICM, which is best resolved with Chandra , we 
will not consider the XMM-Newton data further here. 


4. SPATIALLY RESOLVED SPECTROSCOPY 

Unless otherwise specified, all spectra were fitted in the 
0.6-3.0 keV band using xspec, with an absorbed APEC 
model with the absorption fixed at the Galactic value 
of Nh = 4.37 x 10 20 cm -2 (Kalberla et al. 2005) and 
the abundance allowed to vary. Spectra were grouped 
with a minimum of 40 counts per energy bin. Anders & 
Grevesse (1989) abundance ratios and AtomDB 2.0.2 
were used throughout. We note that using this version 
of AtomDB gives systematically larger temperatures (by 
roughly 0.05-0.1 keV, significantly larger than our typi- 
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cal statistical error) as compared with AtomDB 1.3.1, 
which was used in Rll. This difference arises from sig¬ 
nificant changes in the ionization balance, particularly at 
the low (^0.7 keV) temperatures we consider here (see 
Fig. 3 in Foster et al. 2012). Thus, we expect a system¬ 
atic offset between the spectroscopic quantities that we 
derive here as compared with Rll. 


4.1. Spectroscopic Maps 

To study the thermal structure of the ICM, we gen¬ 
erated a smoothed temperature map using the method 
described in Randall et al. (2008), and also used in Rll. 
For each pixel in the temperature map, spectra were ex¬ 
tracted from a surrounding circular region that contained 
1500 source counts. The temperature map pixel value 
was then set equal to the best fitting temperature from 
fitting to these spectra. Since the extraction regions are 
generally larger than the pixels, nearby pixels are corre¬ 
lated with one another, such that the temperature map 
is effectively smoothed on the scale of the local extrac¬ 
tion region size. Some advantages of this method are 
that there is no a priori assumption about the thermal 
structure of the gas, and experience has shown that such 
maps can resolve real structure on scales that are smaller 
than the size of the extraction regions. 

For comparison, we also constructed temperature maps 
using the contour binning method developed by Sanders 
(2006). Here, extraction regions are defined based on sur¬ 
face brightness contours. Regions were defined to achieve 
a signal-to-noise ratio of 38 (corresponding to ^1,400 
counts in bright regions) in each region. This method as¬ 
sumes that the temperature distribution follows the sur¬ 
face brightness distribution (which is frequently, but not 
always, the case). Some advantages of this method are 
that each extraction region is independent of the others 
(as long as the region widths are large compared with the 
local PSF), and it is relatively cheap computationally (as 
compared with the smoothed temperature map method). 
A comparison of the smoothed and contour binned tem¬ 
perature maps gives an indication of the extraction re¬ 
gion size as a function of position in the smoothed map. 

Smoothed temperature maps over a wider area and a 
higher resolution map of the core are shown beside the 
contour binned maps in Figures [5] and [ 6 | respectively. 
Both the smoothed and contour binned temperature 
maps reveal relatively hot emission associated with the 
bright rims surrounding the inner cavities and just inside 
the prominent 10 kpc surface brightness edges. These 
features were identified as shocks in Rll. The maps 
suggest additional temperature rises associated with the 
outer 30 kpc surface brightness edges, particularly to the 
NW. The te mpera ture profiles across these edges are pre¬ 
sented in § |4(2.2| The temperature maps also show a 
plume of cooler gas, roughly extending along the line de¬ 
fined by the cavities, from SW to NE. In Rll, we showed 
that this plume is consistent with arising from cool gas 
that has been lifted by the buoyantly rising X-ray cavi¬ 
ties. 

We derived pseudo-pressure and pseudo-entropy maps 
from the smoothed temperature maps. The fitted APEC 
normalization for each spectral map pixel is proportional 
to the volume integral of the square of the electron den¬ 
sity along the line of sight. The square root of the nor¬ 
malization therefore gives an average “projected” density 


along the line of sight. We define the pseudo-pressure and 
pseudo-entropy as kTA 1 / 2 and &TA -1 / 3 , respectively, 
where kT is the fitted projected temperature and A is 
the normalization scaled by the area of the extraction 
region. The wide field and core maps are shown in Fig¬ 
ures [7] and [ 8 ] The pseudo-pressure maps reveal promi¬ 
nent pressure jumps across the 10 kpc shock fronts and 
increased pressure in the bright central rims (i.e., the 
1 kpc shock). There is no obvious equivalent structure 
in the pseudo-entropy maps. This is not unexpected for 
weak shocks. For the strongest shock we detect, with 
M = 1.78, the Rankine-Hugoniot shock jump conditions 
for a 7 = 5/3 gas give a pressure jump at the shock front 
of a factor of 3.7, but a jump in the entropy index of only 
a factor of 1.12. Thus, entropy jumps are expected to be 
relatively difficult to detect, particularly in projection. 

The core pseudo-entropy map (Figure | 8 | shows some 
filamentary structure in the center, on the scale of the 
bright shocked rims. However, unlike the features in 
the pseudo-pressure map, these features do not directly 
trace the bright rims. Furthermore, the entropy is lower 
in these filaments as compared with the surrounding re¬ 
gions, in contrast to the small increase in entropy that 
is expected at the shock fronts. These structures are 
likely formed as the central low-entropy gas is pushed 
and pulled away from the center as the inner cavities 
expand and rise buoyantly. 

The central regions of the core temperature and 
pseudo-entropy maps are shown with the Ha contours 
from Rll overlaid in Figure [9j As shown in Rll, the 
Ha emission follows the plume of cool X-ray emitting 
gas up to the location of the intermediate cavities, in¬ 
dicating that these cavities have uplifted cold phase gas 
from the core as they buoyantly rise. The Ha emission 
also clearly traces the SW inner cavity, suggesting that 
the central cold-phase gas is dynamically coupled to the 
cavity and is being pushed out as the cavity expands. 
The SE 1 kpc shock appears to be passing through the 
Ha filament (at least in projection) without destroying 
it. Finally, there is a correlation between the Ha emis¬ 
sion and the detailed, small-scale, filamentary structure 
of the central low-entropy gas. This is consistent with 
the above suggestion that these entropy features likely 
arise from gas that has been displaced from the core by 
the cavities, rather than from entropy changes driven by 
the central shocks. 

Since each extraction region in the above maps con¬ 
tained only 1500 net counts, the abundance values deter¬ 
mined by these fits were not tightly constrained, and the 
corresponding abundance maps show little to no struc¬ 
ture. Abundance measurements in ^0.7 keV gas are chal¬ 
lenging due to the large number of closely spaced emis¬ 
sion lines, and the resulting degeneracy between line and 
continuum emission at Chandra’s energy resolution. We 
found that, for N5813, achieving 1 -cr errors on the order 
of 5-10% for the abundance required roughly 30,000 net 
counts. Fortunately, the new deeper observations pro¬ 
vided well over one million photons appropriate for use 
in spectral fitting, allowing us to map the abundance dis¬ 
tribution in the gas. A smoothed abundance map, with 
30,000 net counts per extraction region, is shown in Fig¬ 
ure [lOl 

The abundance map shows an apparent plume of low 
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abundance gas that corresponds to the plume of up¬ 
lifted cool gas seen in Figure |5] (roughly 40 kpc long, 
extending from SW to NE). It has long been recognized 
that fitting multi-temperature spectra with a single ther¬ 
mal model can give anomalously low abundance values 
(the so-called “Fe-bias” effect, Buote 1999). To test for 
this effect, we fit multi-temperature models to the ap¬ 
parent central minimum and apparent low-metallicity 
NE plume. We find that the plume is well described 
by a two-temperature model, with kT\ ow « 0.40 keV, 
^Thigh ~ 0.75 keV, and an abundance (set equal for each 
temperature component) of Z p i ume ~ 45% solar, consis¬ 
tent with the surrounding ICM. In contrast, we find that, 
while including multiple temperature components for the 
central region (in a 2.2 x 3.8 kpc ellipse oriented NE to 
SW, along the plume) does raise the fitted abundance, 
it is still significantly lower than the abundance of the 
surrounding ICM, at Z center « 37% solar. Thus, we con¬ 
clude that the apparent low abundance in the extended 
plume is a fitting artifact arising from the Fe-bias effect, 
and that the measured abundance in the plume is consis¬ 
tent with that of the surrounding ICM, while the appar¬ 
ent abundance minimum within the central few kpc can 
not be easily explained by this effect. We note that cen¬ 
tral abundance dips that apparently cannot be explained 
solely by projection effects have been seen in other sys¬ 
tems (e.g., Blanton et al. 2003; Panagoulia et al. 2013). 
In addition to the apparent low abundance plume, there 
appear to be increases in the projected abundance at the 
location of the shock fronts, in contradiction with what is 
expected from the Fe-bias effect. A more detailed study 
of these features and the ICM abundance in general in 
N5813 will be presented in an upcoming paper. 

4.2. Radial Profiles 
4.2.1. Azimuthally Averaged Profiles 

We produced projected radial profiles by fitting spectra 
extracted from concentric annuli, centered on the cen¬ 
troid of the diffuse emission at large radii. Each annu¬ 
lus was chosen to be at least 1" wide and contain at 
least 1000 net counts in the 0.6-3.0 keV band. The inner 
annuli were limited by the size constraint (and there¬ 
fore contained more than 1000 net counts), whereas the 
outer annuli were wider and contained roughly 1000 net 
counts. The resulting temperature profile is shown in 
the top panel of Figure [Til The positions of the shock 
fronts (determined in § p| are marked by vertical lines. 
Even in the projected temperature profile in circular an¬ 
nuli (which effectively smooths out the elliptical shock 
fronts), temperature enhancements associated with each 
shock are clearly visible. The central 1.5" radius region, 
which contains the central AGN, has been excluded from 
the fit to the central region. The central temperature rise 
is due to the presence of the inner shock rims surrounding 
the inner cavities. 

As in Rll, we determined the 3D structure of the ICM 
by performing an “onion peeling” deprojection analy¬ 
sis (Fabian et al. 1981) using concentric annuli, each 
of which contained at least 8000 net counts (in projec¬ 
tion) in the 0.6-3.0 keV band. We fit each annulus with 
an absorbed APEC model, with the abundance fixed at 
50% solar (consistent with results from spectral fits to 
the total diffuse emission and with typical values in the 


projected abundance map shown in Figure Jl0|). We as¬ 
sume spherical shells of uniform density. Given the rich 
morphology of N5813, with elliptical shock edges, cavi¬ 
ties, and other features, it is clear that the assumption 
of spherical symmetry does not strictly apply. Thus, the 
deprojected profiles include some level of systematic un¬ 
certainty associated with the assumed geometry. The 
deprojected temperature, density, pressure, and entropy 
profiles are shown in Figure [TTJ We see that the de¬ 
projected temperature is generally consistent with the 
projected temperature. We interpret this an indicating 
that the scale of the azimuthally averaged temperature 
gradient is generally small compared with the sizes of 
our extracted annuli, so that projection effects in any 
given annulus (which are dominated by emission from 
adjacent annuli) are relatively small. The small discon¬ 
tinuous jumps, or “kinks”, in the density, pressure, and 
entropy profiles around 11 kpc are associated with the 
sharp surface brightness edges to the NW and SE from 
the intermediate shock. 


4.2.2. Shock Front Profiles 

To confirm the temperature rises across the surface 
brightness edges seen in Figures [5] and [6j we derived tem¬ 
perature profiles in annular bin sectors across each of the 
edges. The profiles were centered on the center of curva¬ 
ture for each feature (such that there was a bin bound¬ 
ary that traced the surface brightness edge), then the 
distances were corrected to measure the distance from 
the center adopted for the azimuthal profiles presented 
4X1] 


The temperature profiles across the 10 kpc and 30 kpc 
shock fronts are shown in Figures [ 12 ] and [l3j respective! 

)nt positions overls 


with the fitted shock front positions overlaid (see § |5j) 
Clear temperature rises are seen at each front, confirm¬ 
ing that these features are shocks. In the case of the 
10 kpc shocks, the deeper observations allow us to mea¬ 
sure the temperature profiles with greater angular reso¬ 
lution (on the scale of a few hundred pc) as compared 
with Rll. In addition to the overall shift in normaliza¬ 
tion due to the updated AtomDB (see § [4|, the peak 
change in temperature across the fronts is larger as com¬ 
pared with the more coarsely binned profiles in Rll. In 
contrast, the shock Mach numbers that we derive here 
(§ [5]) are consistent with those derived in Rll, since the 
surface brightness profile requires fewer counts per radial 
bin, allowing for smaller bins. This is consistent with our 
suggestion in Rll that the temperature jumps we found 
there were smaller than what was expected based on the 
Mach numbers derived from the surface brightness jumps 
due to projection effects. 

In some cases, there appears to be shock heated gas 
just outside the positions of the shock front (e.g., the two 
bins to the right of the edge in Figure [l 2 j left). There 
are at least two factors that can contribute to this effect. 
First, there is a small uncertainty in the fitted locations 
of the shock fronts. Second, we show in § [5] that the 
shock fronts have finite widths (on the order of 0.5 kpc 
for the 10 kpc shock). When fitted with a discontinuous 
model, the preferred position of the edge will naturally 
be at the center of the smoothed front. 


5. STRUCTURE OF THE SHOCK FRONTS 
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To quantitatively characterize the structure of the 
shock fronts, we applied the edge fitting method de¬ 
scribed in Rll. We extracted the 0.3-3.0 keV surface 
brightness profiles across each of the three shocks (at 
roughly 1 kpc, 10 kpc, and 30 kpc) to the SE and the 
NW, where the edges are sharpest. We excluded the 
NW side of the 1 kpc shock as this region shows a compli¬ 
cated morphology, with the shocks driven by the NE and 
SW inner cavities just beginning to overlap (Figure [lj, 
and is thus not well described by our simple spherically 
symmetric model. Although the 1 kpc SE shock edge 
appears relatively regular, we cannot rule out the possi¬ 
bility of some systematic bias in our results due to over¬ 
lapping cavity rims. The level of this potential bias is 
difficult to characterize without detailed numerical sim¬ 
ulations, although we do not expect it to strongly affect 
our results. Each surface brightness profile was centered 
on the center of curvature for the corresponding edge, 
which was not coincident with the location of the AGN 
or with th e adop ted center for the azimuthal profiles de¬ 
rived in § |4.2.1| These surface brightness profiles were 
converted into integrated emissivity profiles (IEM) using 
the best fitting projected temperature in each bin (deter¬ 
mined from fitting a single absorbed APEC model) with 
the abundance fixed at 50% solar. The IEM profiles were 
then fit by projecting a spherically symmetric discontin¬ 
uous power-law density model. The free parameters of 
this model are the normalization, the inner density slope, 
the outer density slope, and the amplitude and location 
of the density jump. Finally, the density jumps were 
converted into shock Mach numbers using the standard 
Rankine-Hugoniot shock jump conditions for a 7 = 5/3 
gas. The fit to the NW 10 kpc shock is shown in Fig¬ 
ure [Ml The results of all of the fits are summarized in 
Table 2] (note that all radii are scaled to give the dis¬ 
tance from the cente r adopted for the azimuthal profiles 
presented in § |4.2.1[ rather than the distance from the 
center of curvature of each front). We find Mach num¬ 
bers of roughly 1.8, 1.5, and 1.2 for the 1 kpc, 10 kpc, 
and 30 kpc shocks, respectively. These results confirm 
the presence of the outer shock at ~ 30 kpc, only hinted 
at in the shallower observations analyzed in Rll. This 
shock is associated with the same AGN outburst that 
inflated the outermost cavity pair ^ 50 Myr ago. 

The IEM profiles of both the NW and SE 10 kpc shock 
fronts show edges that are less sharp than the discontin¬ 
uous power-law model. We therefore fit each edge with 
the same model but smoothed with a Gaussian, with the 
width of the Gaussian (a) as an additional free parameter 
of the fit. The fit for the NW 10 kpc shock front is shown 
along with the discontinuous model fit in Figure [Mj The 
smoothed model is a significantly better fit to the data 
for both the SE and NW 10 kpc shock fronts, with both 
probabilities < 0.04% based on the F test. For the SE 
1 kpc shock, the smoothed model is a marginally better 
fit, with a probability of ^ 9%, while for the 30 kpc shock 
both models describe the data equally as well. The 1 kpc 
SE and both 10 kpc shock fronts are shown to have fi¬ 
nite widths at greater than 90% confidence. Only upper 
limits can be placed on the widths of the 30 kpc shock 
fronts, where the surface brightness is relatively low. 

The shock front edges are expected to be blurred by 
Chandra’s PSF, which varies across the field. Table [ 2 ] 
gives the size of the 50% encircled energy fraction (EEF) 


region at the location of each shock, assuming a 0.74 keV 
plasma with an abundance of 50% solar and Galactic ab¬ 
sorption, and the fitted edge widths. The 90% EEF is 
about twice as wide as the 50% EEF for this model. For 
the SE 1 kpc shock front, the edge width is on the or¬ 
der of the size of the PSF. Thus, the measured width of 
this edge may not be intrinsic to the front. However, for 
the NW and SE 10 kpc shock fronts, the 90% confidence 
lower limit on a is larger than the local PSF. The differ¬ 
ence is particularly significant for the NW 10 kpc shock, 
where the PSF is smaller (the optical axis is located NW 
of the central core). Furthermore, we note that the best 
fit shock width is somewhat larger (consistent within the 
errors) for the NW 10 kpc shock front, as compared with 
the SE 10 kpc front, despite the fact that the PSF is 
more than twice as large in this region. This shows that 
the measured widths are not driven by PSF blurring. 

We conclude that, in the case of the 10 kpc shock, the 
edge widths given in Table [2] are likely intrinsic to the 
edges themselves and not instrumental artifacts. In the 
case of the 30 kpc shock, only upper limits on the shock 
widths can be placed. The intrinsic widths may simply 
be due to irregularities in the shape of the shock front, 
due to variations in the shock strength with azimuthal 
angle created as the shock propagates through a non- 
uniform and/or turbulent ICM. To test this effect, we 
fit the surface brightness profiles, as above, in sectors of 
varying angular width for the SE and NW 10 kpc shock 
fronts, excluding very wide sectors that include the X- 
ray cavities. We found consistent density jumps for every 
case that we tried, although for very narrow sectors the 
uncertainties were large. Thus, while we cannot rule out 
the possibility that the finite shock front widths that we 
find are due to deformations in the shock fronts (possibly 
at smaller angular scales than we are able to measure), 
we do find that determining the shock widths in more 
narrow sectors does not result in significantly smaller 
widths. Possible implications of finite shock widths and 
the impact of a turbulent ICM are discussed further in 

§ [673 

6. DISCUSSION 

6 . 1 . Heating, Cooling, and Feedback in the ICM 

Shocks are expected to heat the ICM as they propa¬ 
gate. While there is a relatively large temperature jump 
at the shock front, this increase is transient as the gas 
subsequently expands and cools. The relevant quantity 
for lasting heating of the gas is the change in entropy. 
The equivalent amount of heat energy A Q imparted to 
the gas by a shock due to a change in entropy AS is given 
by 

AQ ~ TAS = EA In ( 1 ) 

where E = Cy T is the total thermal energy of the gas, 
7 is the adiabatic index (taken to be 5/3), and p and 
p are the pressure and density, respectively. Thus, each 
shock contributes a fraction A In ^ of the total thermal 
energy in the gas. In the case of weak shocks, this frac¬ 
tion is small (between 0.4-12% for the shocks we detect 
in N5813). To compare with radiative cooling, we are 
interested in the cumulative effect of shock heating per 
local cooling time in the gas (i.e., the number of shocks 
per local cooling time and the total amount of heat en- 
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ergy deposited in the gas). We approximate the local 
cooling time just outside each shock edge as the time it 
would take the gas to radiate away all of its thermal en¬ 
ergy based on its deprojected emissivity. Based on the 
shock ages given in Table [2j and keeping in mind that 
these ages are expected to be underestimates for young 
shocks and overestimates for older shocks, we assume an 
outburst repetition rate of one every 20 Myr. 

The results from the above comparison of the shock 
heating and radiative cooling rates are summarized in Ta¬ 
ble [3] Given the various statistical and systematic uncer¬ 
tainties associated with, e.g., the shock Mach numbers, 
the shock ages and outburst repetition rate, the depro¬ 
jected cooling luminosity, etc., these numbers should be 
considered to be rough estimates. Nevertheless, there is 
good agreement between the number of shocks expected 
per cooling time (column 3) and the number required 
to completely offset radiative cooling (column 4). We 
conclude that shock heating alone is sufficient to balance 
radiative cooling in the gas out to at least ~30 kpc. This 
extends our earlier result in Rll, where we found a bal¬ 
ance between heating and cooling out to ~10 kpc. At 
30 kpc, the cooling time (^ 2 Gyr) is approaching the 
expected age of the system, indicating that there is no 
need for additional heating beyond this radius to com¬ 
pletely balance cooling and explain the lack of a strong 
cooling flow in this system. 

We can estimate the average cavity power P cav by di¬ 
viding the total enthalpy of the cavities (estimated as 
4 pV for each cavity, for a total of 5.7 x 10 57 erg s -1 ) 
by the age of the oldest cavity pair (^ 10 8 yr; see 
Rll, Table 3). This gives P cav w 1.8 x 10 42 erg s _1 . 
Spectral fitting gives a total X-ray luminosity of Lx ~ 
5.5 x 10 41 erg s -1 within a radius enclosing all of the 
cavities (~ 26 kpc). Thus, in N5813 we find that a few 
times the pV work required to inflate the cavities is in 
principle large enough to offset radiative cooling, as has 
generally been found in other systems (e.g., Birzan et al. 
2004; Rafferty et al. 2006; Nulsen et al. 2007; Hlavacek- 
Larrondo et al. 2012). 

Since we have shown that, in N5813, shock heating 
alone is sufficient to balance cooling, one might ask 
whether the total heating from cavities and shocks leads 
to an “overheating” of the ICM, breaking the balance 
between heating and cooling that is required for AGN 
feedback to operate over long time scales. Although cav¬ 
ities are expected to deposit some of their enthalpy in 
the ICM as they rise buoyantly (McNamara & Nulsen 
2007), this heating is not directly observed, and the de¬ 
tails on how and where cavities heat the ICM are poorly 
understood. Furthermore, cavities can only directly heat 
the ICM locally, as they rise radially, whereas roughly 
isotropic heating must take place close to the central 
AGN to regulate feedback. Mixing may help to distribute 
locally heated gas (e.g., Zhuravleva et al. 2014), although 
the gas motions must be small enough to not completely 
destroy the observed X-ray cavities. Shocks, in contrast, 
heat the gas roughly isotropically and most strongly at 
smaller radii, and the thermal effect of shocks on the ICM 
is, in some cases, directly observed. We suggest that 
AGN outburst shocks may generally play a significant 
role in AGN feedback, particularly at small radii where 
the Mach numbers are higher. The cavities are then free 


to rise buoyantly, and release their energy to heat the 
ICM at larger radii. The fact that cavities are observed at 
larger radii demonstrates that they do not necessarily re¬ 
lease their enthalpy to the ICM at small radii, consistent 
with this picture. Weak shocks are generally difficult to 
detect, since the shock fronts are thin and easily masked 
by projection effects. They are more easily detected in 
N5813 due to its proximity, low gas temperature, and 
regular morphology. The fact that one of the handful 
of examples of confirmed outburst shocks known in clus¬ 
ters, in Abell 2052 (Blanton et al. 2011), was confirmed 
by a measured temperature jump only after very deep 
Chandra observations and careful deprojection is consis¬ 
tent with this scenario. Similarly, no temperature jump 
is found associated with the presumed outburst shock 
in the Perseus cluster using extremely deep Chandra ob¬ 
servations, even though it is the X-ray brightest cluster 
(Graham et al. 2008). 

6.2. AGN Outburst History 

Measurements of the cavities and shocks allow us to 
estimate the total mechanical outburst energy for each 
of the three outbursts in N5813. Since the mechanical 
luminosity dominates the total AGN energy output in 
N5813 (Rll), as is generally the case for kinetic mode 
AGN feedback (Fabian 2012), this gives us information 
on the total outburst energy history of the AGN. We take 
the cavity internal energy to be 3pF, where the pressure 
p is taken from the azimuthally averaged pressure pro¬ 
file. The assumed cavity dimensions and derived internal 
energies are given in Table [4] We use the cavity inter¬ 
nal energy rather than the enthalpy (which is the sum of 
the internal energy and the work required to inflate the 
cavity, and roughly equal to 4 pV for cavities filled with 
relativistic particles) since the work done during cavity 
expansion goes into driving shocks, the energy of which 
we account for separately. As in Rll, we estimate the 
outburst energy in shocks as 

E s Mp 1 V s (p 2 /pi~l), (2) 

where pi and P 2 are the pre- and post-shock pressures, 
respectively, and V s is the total volume enclosed by each 
roughly ellipsoidal shock surface. Although equation [2] 
is expected to slightly underestimate the shock energy 
as the Mach numbers were larger at earlier times (giving 
larger pressure jumps), in Rll we showed that this esti¬ 
mate agrees well with shock energies derived from simple 
ID hydrodynamic simulations of a central point explo¬ 
sion by matching the simulated surface brightness pro¬ 
files with the observations. The total mechanical energy 
output of each outburst is simply the sum of the total 
cavity and shock energies. We note that the ratio of the 
energy in cavities to that in shocks is between 0.15-0.3 
for each outburst, and that in the case of the central 
outburst, where the cavities occupy a significant fraction 
of the volume enclosed by the 1 kpc shock surface, this 
ratio is somewhat lower than the expected value of ~ 1. 
This is likely a result of underestimating the central gas 
pressure, and hence the central cavity energies, due to 
the complicated morphology in this region, which leads 
to systematic uncertainties associated with projection ef¬ 
fects and a lack of knowledge on the volume filling factor 
of the X-ray emitting gas. However, we note that such 
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an underestimate would only moderately affect our total 
outburst energy, well within the factor of a few at which 
we expect our simple estimates to be accurate, and hence 
would not affect our main conclusions. 

Using the above estimates, we find total outburst en¬ 
ergies of 8.5 x 10 56 erg, 9.9 x 10 57 erg, and 8.9 x 10 57 erg 
for the outbursts associated with the 1 kpc, 10 kpc, and 
30 kpc shocks, respectively. As in Rll, we find the out¬ 
burst energy of the 1 kpc outburst to be more than a fac¬ 
tor of ten less than that of the 10 kpc outburst. Based on 
this, we concluded that either the 1 kpc outburst is ongo¬ 
ing, with the AGN actively inflating the cavities, or that 
the mean jet power varies significantly over timescales 
on the order of 10 Myr. Here, we find remarkable agree¬ 
ment between the total outburst energies of the 10 kpc 
and 30 kpc outbursts. This suggests that the mean out¬ 
burst (or jet) power, averaged over the outburst interval 
of a few 10 Myr, is fairly stable, and that the total en¬ 
ergy of the 1 kpc outburst is lower because this outburst 
is ongoing. 

In Rll, we estimated the outburst power as the to¬ 
tal energy divided by the age of the outburst. Here, 
we instead estimate this power as the total energy di¬ 
vided by the duration of each outburst, which we take 
to be the shock age for the youngest ou tbur st (Table [2| 
and the outburst interval (^ 20 Myr, § 6.1| ) for the re¬ 
maining outbursts. This gives mean outburst powers of 
1.6xl0 43 erg s -1 , 1.6xl0 43 erg s -1 , and 1.4xl0 43 erg s -1 
for the 1 kpc, 10 kpc, and 30 kpc outbursts, respectively. 
Thus, we find the mean outburst power is indeed roughly 
constant, certainly within a factor of a few, which is the 
level at which we expect our rough estimates to be accu¬ 
rate. We note that in Rll we found a lower power for the 
1 kpc outburst. This difference arises mainly from using 
results from our numerical model in Rll, which predicts 
a lower shock energy and larger shock age as compared 
with the estimates based on observations alone (see Ta¬ 
ble 2 in Rll). The source of this difference is likely to be 
a systematic error that arises when applying the spheri¬ 
cally symmetric model to the 1 kpc shock, which clearly 
has a more complicated morphology (see Figures [l]&[6|. 

We stress that our result on the consistency of tne out¬ 
burst power only applies to the mean power averaged 
over the outburst interval of a few 10 Myr. AGN lu¬ 
minosities are known to vary by up to several orders of 
magnitude over very short (observable) timescales (e.g., 
Harris et al. 2009). However, in principle the outburst 
signatures in N5813 could be created with a constant 
jet power. In this case, the cavities expand rapidly and 
drive shocks just after they are formed. The expansion 
rate drops as the fractional energy input rate decreases, 
and the shocks separate from the cavities and propagate 
outwards. As the cavities rise buoyantly, they eventually 
disconnect from the central jets (once the buoyant speed 
exceeds the expansion rate), which then begin inflating 
a new pair of cavities and the process repeats (similar 
models are discussed, e.g., in Fabian et al. 2003). There¬ 
fore, we use the term “outburst” to refer to the creation 
of a cavity pair and its associated shock, rather than a 
rapid increase in jet power, since these features are con¬ 
sistent with a kinetic jet power that either also cycles in 
outbursts or is roughly constant. 

We conclude that the most likely scenario is that N5813 
is in a “steady state” kinetic feedback mode, with the 


AGN outbursts roughly equally spaced in time, and each 
outburst depositing a similar amount of total energy into 
the ICM, such that the mean AGN outburst power is 
roughly constant. The lower total power of the 1 kpc out¬ 
burst is an indication that this outburst is ongoing, and 
has yet to deposit the bulk of its energy in the ICM. The 
detection of 1.4 GHz radio emission filling the inner cav¬ 
ities (from young, high-energy, non-thermal particles in¬ 
jected by the AGN), along with shock heated cavity rims 
that have not yet separated from the cavities (as they 
have clearly done for the intermediate and outer cavi¬ 
ties), are consistent with this interpretation (Rll). This 
is also consistent with results from Allen et al. (2006), 
who find evidence for accretion flows around central AGN 
that are stable over a few million years. The fact that 
N5813 shows little evidence of a recent merger or other 
complicated ICM “weather”, probably contributes to it 
being able to maintain such steady state feedback over 
long timescales. 

6.3. Transport Processes in the ICM 

Fits to the surface brightness profiles across the shock 
fronts reveal that the SE and NW 10 kpc fronts have 
non-zero widths that cannot be explained by Chandra’s 
PSF (§[§. Measurements of the widths of surface bright¬ 
ness edges in cluster X-ray observations have been used 
to place constraints on transport processes in the ICM by 
comparing these widths with the local collisional mean 
free paths of the particles. Edges that are significantly 
narrower than the particle mean free path indicate that 
Coulomb diffusion is suppressed across the edge. In the 
case of cold fronts in clusters, the front widths are found 
to be significantly smaller than the particle mean free 
paths (e.g., Vikhlinin et al. 2001a; Russell et al. 2012), 
implying that Coulomb diffusion is suppressed across 
the fronts. However, in cold fronts the suppression of 
transport processes (and Kelvin-Helmholtz instabilities) 
is likely due to magnetic draping, where magnetic field 
lines are stretched along the cold front edge as the rela¬ 
tively cool, dense gas moves through the ICM (Vikhlinin 
et al. 2001b; Asai et al. 2004, 2005, 2007; Lyutikov 2006; 
Dursi & Pfrommer 2008, ZuHone et al. 2011). In this 
case we consider shock fronts which propagate through 
the ICM, and therefore magnetic draping is not expected 
to occur. 

In the case of merger shocks, Markevitch & Vikhlinin 
(2007) find a bow shock width of roughly 35 kpc, on the 
order of the local particle mean free path, for IE 0657- 
56 (the Bullet cluster), although this width is only 
marginally preferred over a zero width front. Russell 
et al. (2012) examine shock widths in the merging clus¬ 
ter Abell 2146, which contains a leading bow shock 
and a trailing reverse shock. They find finite width 
shocks in each case, with the reverse shock width be¬ 
ing significantly smaller, and the bow shock width be¬ 
ing marginally smaller, than the particle mean free path. 
They conclude that transport processes are suppressed 
across these merger shocks. In the case of AGN outburst 
shocks, Croston et al. (2009) consider the thickness of 
the northeastern shock in Centaurus A, but their results 
are inconclusive. 

It is therefore of interest to compare our measured 
shock widths with the Coulomb mean free path of parti¬ 
cles in the ICM. In the region of the shock fronts there 
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are four relevant mean free paths: that in the pre-shock 
region, Ai n , that in the post-shock region, A out , that of 
particles crossing the front from the post-shock to the 
pre-shock region, Ai n -K> u t> and that of particles crossing 
the front from the pre-shock region to the post-shock re¬ 
gion, Aout^in- In our case, the largest and most relevant 
mean free path is Ai n _^ 0 ut ? which is given by 



( 3 ) 


where n e , out is the pre-shock electron density, x = 
Tin/T 0 ut, Ti n is the post-shock gas temperature, Tout 
is the pre-shock gas temperature, G(y) = [4>{y) — 
ycf)'(y)\/2y 2 , and cj)(y) is the error function (Spitzer 1962). 
When calculating Ai n _^ out , pre-shock temperatures were 
taken from the projected temperature profiles just out¬ 
side the shock fronts. Post-shock temperatures were cal¬ 
culated from the pre-shock temperatures by applying the 
Rankine-Hugoniot jump conditions, rather than taken 
from the observed projected post-shock temperatures. 
We note that this is a conservative assumption from the 
perspective of finding a mean free path that is smaller 
than the shock width, as it gives a higher post-shock tem¬ 
perature (that is not diminished by projection effects), 
which gives a larger value for Ai n ^ 0 ut- 

The relevant particle mean free path is given for each 
shock front in Table [2j The uncertainty in the mean free 
path is dominated by the uncertainty in the pre-shock 
temperature, which is on the order of 10%, leading to a 
mean free path uncertainty of roughly 20%. In all cases, 
the particle mean free path is significantly smaller than 
the shock width, on the order of ten times smaller than 
the 90% confidence lower limit for the 1 kpc and 10 kpc 
shocks (where lower limits on the shock widths can be 
placed). For a weak shock, the shock width is expected to 
be roughly w « A/(M — 1), where A is the effective parti¬ 
cle mean free path (Landau & Lifshitz 1987; McNamara 
& Nulsen 2007), or 1.4-2 times A for the weak shocks 
we consider here. Thus, the apparent shock widths are 
too large to be explained by particle diffusion alone. We 
note that this result likely precludes using these shock 
fronts to place accurate constraints on the electron-ion 
equilibration timescale, as has been done, e.g., by Russell 
et al. (2012). 

As noted in § [5j a measured finite shock width may 
arise due to deformations in and broadening of a shock 
front as it propagates through a “clumpy” and/or turbu¬ 
lent ICM. Nulsen et al. (2013) provide an estimate of the 
expected shock width due to turbulence as a function of 
radius as the shock propagates through a uniformly tur¬ 
bulent ICM. The shock width (w, defined as the rms of 
the displacement of sections of the shock front due to tur¬ 
bulence) is determined by the distance traveled (r s ), the 
shock speed u s , the coherence length of the turbulence 
(£), and the rms turbulent speed (cr t ). Here, we can in¬ 
vert this relation to find the turbulent speed implied by 
the observed shock width: 


at 


WVs_ 

VrJ-' 


( 4 ) 


The rms turbulent speed implied by measurements of 
each shock front we consider are given in Table [2j Fol¬ 
lowing Nulsen et al. (2013), we take the coherence length 


to be i ~ O.lr (Rebusco et al. 2005). Remarkably, 
the implied turbulent speed is consistent with roughly 
70 km s -1 in each case. We note that, in principle, 
the turbulent speeds calculated here are upper limits on 
the true ICM turbulence, as this calculation does not 
include the contribution of density inhomogeneities to 
shock broadening, projection effects, and the effects of an 
irregular shock front geometry (which is approximated as 
spherically symmetric in sectors), all of which will act to 
increase the apparent shock width. However, as the cal¬ 
culation is only approximate, our results should not be 
taken as hard upper limits. 

The implied turbulent velocity of roughly 70 km s -1 is 
reasonable compared with results from simulations (100- 
300 km s -1 ; Lau et al. 2009) and observations (e.g., de 
Plaa et al. 2012; Sanders & Fabian 2013; Zhuravleva 
et al. 2014). While this value is on the low end of the 
reported range, we note that N5813 is a relatively iso¬ 
lated galaxy group (Rll), which has a regular morphol¬ 
ogy and shows no signs of a recent major merger. Thus, 
the level of merger-driven turbulence is expected to be 
low. Additionally, we note that observational measure¬ 
ments will naturally be biased towards detecting systems 
with larger turbulent velocities. In several cases only up¬ 
per limits can be placed, with all reported upper limits 
significantly larger (a few to several 100 km s -1 ) than 
the ~ 70 km s -1 detection we report here (Sanders et al. 
2011; Bulbul et al. 2012b; Sanders & Fabian 2013). Fi¬ 
nally, this value is consistent with the range of turbulent 
velocities of 43-140 km s -1 recently reported by Zhu¬ 
ravleva et al. (2014) for the (relatively low mass) Virgo 
cluster. 

By considering the effects of resonant scattering 
on Fe XVII line ratios, de Plaa et al. (2012) con¬ 
strain the turbulent velocity in N5813 to be 140 < 
kturb < 540km s -1 , significantly larger than our value 
of 70 km s -1 . However, we do not view this as a seri¬ 
ous discrepancy for several reasons. First, as mentioned 
above, the limits provided by equation [4] are expected to 
be rough estimates. In particular, the coherence length 
i is unknown and the only hard constraint is that i < r. 
Second, as noted by the authors themselves, the atomic 
data used to calculate the line ratios in de Plaa et al. 
(2012) suffer from significant systematic uncertainty. For 
example, by comparing results derived with the SPE^xpl 
code versus ATOMDBpq they find inferred turbulent ve¬ 
locities for N5813 that differ by almost a factor of two. 
Finally, we note that there are plausible physical origins 
for the higher level of turbulence found by de Plaa et al. 
(2012). For example, they place constraints for the cen¬ 
tral r < 5 kpc only, where the turbulent velocity may be 
larger than the total ICM average due to the influence 
of the AGN, whereas the estimates given here are aver¬ 
ages over the entire region interior to each shock. The 
trend of the inferred turbulent speed decreasing with ra¬ 
dius (Table [2h , although not statistically significant, is 
consistent with this picture. It is interesting to note that 
the turbulent velocity inferred from the width of the cen¬ 
tral r ~ 1 kpc shock is consistent with the range given 
by de Plaa et al. (2012), although, as noted above, the 
errors are large and the effects of the larger relative PSF 

11 http://www.sron.nl/spex 

12 http://www.atomdb.org 
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are unclear. Additionally, we note that the resonant line 
ratio may also be affected by bulk flows in the ICM (e.g., 
the radial gas flows behind the shock fronts), possibly 
leading to an over-estimate of the turbulence. 

We conclude that, while the measured shock widths 
are too large to be explained by Coulomb diffusion, they 
are consistent with arising from the shocks propagating 
through a turbulent ICM. The implied ratio of turbu¬ 
lent to thermal pressure support, estimated as the ratio 
of the turbulent energy to the thermal energy (de Plaa 
et al. 2012), is on the order of a few percent (< 5%) 
for a turbulent speed of 70 km s -1 . For each shock we 
consider, the local thermal diffusion coefficient v ~ Ac s , 
where c s is the sound speed, is significantly smaller (by 
up to an order of magnitude) than the implied turbulent 
diffusion coefficient r] ~ £cr t . Thus, turbulent diffusion 
dominates the microscopic transport in the ICM (at least 
in this case, and likely in other similar-size groups, which 
are expected to have similar bulk ICM properties). 

6.4. Contents of the X-ray Cavities 

Although the X-ray cavities must be close to pressure 
equilibrium with their surroundings, radio observations 
show that generally the pressure support from relativistic 
non-thermal electrons, under the assumption of equipar- 
tition, is insufficient to balance the thermal pressure in 
the ICM (Birzan et al. 2008). Additional pressure sup¬ 
port might be provided by heavy non-thermal particles, 
or by high temperature thermal gas, with deep X-ray ob¬ 
servations providing constraints on the contribution from 
the latter (e.g., Schmidt et al. 2002; Blanton et al. 2003; 
Sanders & Fabian (2007); Rafferty et al. (2013)). Here 
we place similar limits for the low mass galaxy group 
N5813. 

To place lower limits on the temperature of thermal 
gas in the cavities we use a method similar to that em¬ 
ployed in previous studies of other systems (e.g., Sanders 
& Fabian 2007; Rafferty et al. 2013). Due to the complex 
thermal structure of the gas, we do not use comparison 
regions outside of the cavities as in Sanders & Fabian 
(2007). Rather, we place limits based on fits to spec¬ 
tra extracted from regions at the cavity centers, where 
the contrast is highest and the relative contribution from 
high temperature gas inside the cavities is expected to be 
the largest. 

Each spectrum is fit with two absorbed APEC model 
components. The first component models the contribu¬ 
tion from the total projected emission, while the second 
models hotter thermal emission from within the cavities. 
To derive conservative estimates, the temperature, abun¬ 
dance, and normalization of the first component are al¬ 
lowed to vary. We find that within the cavity regions 
this single thermal component is sufficient to adequately 
model the total projected emission (additional thermal 
components do not significantly improve the fits and are 
not well constrained). For the second component, the 
temperature is fixed at some value, the abundance is tied 
to that of the first component, and the normalization is 
allowed to vary. For an assumed cavity volume, the up¬ 
per limit on the normalization gives an upper limit on the 
density. If we further assume that the cavity is in pres¬ 
sure equilibrium with its surroundings, that the pressure 
is equal to the average azimuthal pressure at the cavity 
radius (Figure pT]), and that the thermal gas dominates 


the pressure support inside the cavities, an upper limit 
on the volume filling fraction of the gas in the cavities 
(for the given temperature) is obtained. A lower limit on 
the gas temperature is obtained by increasing the fixed 
temperature of the hot gas until the upper limit on the 
volume filling fraction is equal to unity. 

We concentrate on the cavities with the largest con¬ 
trast and the most regular morphology. Since they are 
closely spaced (and therefore in a similar environment) 
we fit the inner cavities together, assuming a cylindri¬ 
cal geometry with a radius of r = 0.408 (0.488) kpc and 
a depth of d = 0.689 (0.815) kpc for the northeastern 
(southwestern) central cavity. For the intermediate cav¬ 
ities, only the southwestern cavity is considered (with 
r = 1.82 kpc, d = 3.01 kpc) as it has a larger contrast 
and more regular shape compared with the northeastern 
intermediate cavity, which may be split into two cavities 
and may be connected to the northeastern central cav¬ 
ity (Rll). Similarly, we choose the northeastern outer 
cavity over its southwestern counterpart due to its larger 
contrast and more regular shape, and assume an oblate 
spheroidal geometry with major axis a = 5.89 kpc and 
minor axis b = 1.90 kpc in the plane of the sky. We 
find limits on the temperature of any volume filling ther¬ 
mal gas in the cavities of > 3.4 keV, > 16 keV, and 
> 4.5 keV for the inner, intermediate, and outer cav¬ 
ities, respectively (for comparison, the limit placed by 
Rafferty et al. 2013 from deep Chandra observations of 
the galaxy group HCG 62 is 4.3 keV). Our most stringent 
constraint of kT > 16 keV for the intermediate cavities 
implies a gas density in the cavities that is at least 20 
times lower than the azimuthal average at that radius 
(n e < 9 x 10 -4 cm -3 ). 

We note that this method relies on several assump¬ 
tions. For instance, the hot temperature component is 
assumed to arise solely from gas contained within the 
cavities (it may be the case that hot thermal gas sup¬ 
plies only part of the pressure support, while, e.g., heavy 
non-thermal particles or magnetic fields make up the dif¬ 
ference). Additionally, the azimuthally averaged pres¬ 
sure profile, which is used to estimate the mean external 
pressure at each cavity, i s deriv ed under the assumption 
of spherical symmetry (§ |42l| ), which is clearly violated 
at some level for this system (Figure |2|. Furthermore, 
some assumptions about the cavity geometries (or, more 
specifically, their volumes) must be made to calculate 
densities and pressures. The largest source of error is 
likely from the latter assumption. We find that an as¬ 
sumed uncertainty in cavity volume of 20% leads to an 
uncertainty of ^ 10% on the temperature limit (a similar 
result was found by Rafferty et al. 2013). 

6.5. Nature of the Northern Channel 

The residual image shown in Figure [2] reveals a “chan¬ 
nel” of decreased surface brightness to the north, appar¬ 
ently connected to the NE outer cavity. To test the sig¬ 
nificance of this feature, we extracted the surface bright¬ 
ness profile in azimuthal bins roughly 9° wide and 64" 
long across the channel. The profile, shown in Figure p~5| 
shows a clearly significant dip across the channel, indi¬ 
cating that this is a real feature. 

Since the northern channel appears to connect to the 
NE outer cavity, one interesting possibility is that this 
channel is the result of energy (presumably mostly in 
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the form of non-thermal particles) leaking through the 
bright cavity rim and heating the ICM. Indeed, the tem¬ 
perature maps shown in Figure [5] show relatively high 
temperatures in the region of the northern channel, al¬ 
though the correlation between the high temperature re¬ 
gion in the temperature map and the channel in the resid¬ 
ual map is not exact. Alth ough cavities are expected 
to heat the ICM (see § 6.1), such a direct detection of 
this heating has not yet been observed. We therefore 
extracted temperatures both in (wider) azimuthal bins 
across the northern channel, and in identical elliptical re¬ 
gions within and to either side of the channel (at roughly 
the same distance from the center). In neither case did 
we find a significant temperature enhancement associ¬ 
ated with the channel, with temperatures of all regions 
agreeing within roughly 1 -a. We conclude that, although 
it is possible that the northern channel is a result of en¬ 
ergy leaking from the NE cavity and heating the ICM, 
causing the gas in this region to expand to pressure equi¬ 
librium (and subsequently cool adiabatically) and hence 
the density (and surface brightness) to drop, we find lit¬ 
tle evidence to support this interpretation. In particu¬ 
lar, we cannot rule out the possibility that this feature 
is simply due to large scale structure in the ICM, e.g., 
associated with previous mergers or bulk motion of the 
group, even in this relatively relaxed system. Another 
possibility is that this deficit is a left over cavity from 
an even older, fourth outburst that lies off of the axis 
defined by the other three cavity pairs, although the fact 
that it is at least partially within the 30 kpc shock front 
is difficult to reconcile with this interpretation. Deep, 
low-frequency radio observations would be useful to both 
detect the putative low-energy non-thermal particles in 
the NE outer cavity and to see if this emission extends 
out along the northern channel, consistent with it arising 
from a leaking cavity. 


7. SUMMARY 

We have presented results from our analysis of a very 
deep Chandra observation of the galaxy group N5813. 
This observation confirms an outer cavity pair and as¬ 
sociated shock (only hinted at in Rll), giving a total of 
three pairs of roughly collinear cavities, with each pair 
associated with an elliptical shock front with a measured 
temperature jump. The derived Mach numbers are 1.8, 
1.5, and 1.2 for the 1 kpc, 10 kpc, and 30 kpc shocks, 
respectively. These features are cleanly separated signa¬ 
tures from three distinct outbursts of the central AGN. 
We compare the shock heating rate to the radiative cool¬ 
ing rate locally at each shock front, and find that the 
rates closely balance one another. This demonstrates 
that shock heating alone is sufficient to offset gas cool¬ 
ing and enable AGN feedback to operate within at least 
the central 30 kpc (compared with 10 kpc in Rll). We 
suggest that shock heating is likely important at small 


radii in other cool core clusters and groups, but that in 
most cases the shocks are more difficult to detect due 
to higher ICM temperatures, larger distances, and more 
complicated ICM structure, or “weather”, instigated by 
mergers. We find that the total outburst energies of the 
old and intermediate age outbursts are roughly equal, 
while the energy of the young outburst is more than a 
factor of 10 less. However, the mean powers of all of the 
outbursts are roughly the same, within a factor of two, 
indicating that the mean kinetic luminosity of the AGN 
has remained stable for at least 50 Myr. We suggest that 
the young outburst is ongoing, having deposited only a 
fraction of its total energy into the ICM. The proximity 
of the central cavities to the AGN, the relatively high- 
frequency radio emission that fills the cavities, and the 
shock heated rims surrounding the cavities (that have 
not yet separated from the cavities, as is the case for the 
older outbursts) are all consistent with this scenario. We 
find that the 10 kpc (and possibly the 1 kpc) elliptical 
shock front is broadened by ~ 0.4 kpc (~ 0.1 kpc), more 
than ten times the particle mean free path and thus too 
broad to be caused by particle diffusion. While we can¬ 
not rule out broadening due to a clumpy ICM and/or 
projection effects, using rough estimates we find that the 
measured shock widths (which are upper limits for the 
30 kpc shock) are all consistent with broadening due to 
propagation through a turbulent ICM with a turbulent 
speed of ^ 70 km s -1 . This is within, but on the low 
end of, the range of turbulent speeds expected based on 
simulations and other observations, and thus provides a 
robust upper limit on the turbulence since other factors 
may contribute to the total shock broadening. This sug¬ 
gests that transport due to turbulence dominates that of 
particle diffusion throughout the ICM in N5813. Using 
spectral fits in the cavity regions, we place lower limits 
on the temperature of any volume filling gas that could 
completely balance the pressure within the cavity with 
that in the external ICM. Our most stringent limit of 
kT >16 keV comes from the intermediate cavity pair. 
Finally, we find a channel of decreased surface brightness 
extending north from the outer NE cavity. We suggest 
that this feature may be due to energy leaking from the 
cavity and heating the ICM, although we find no tem¬ 
perature enhancement in the channel to support this sce¬ 
nario. 
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TABLE 1 

Chandra X-ray Observations 


Obs ID 

Date Obs 

CCDs Used 

Cleaned Exposure 
(ks) 

5907 

2005 Apr 2 

SI, 

S 2 , 

S3, 

13 

48.1 

9517 

2008 Jun 5 

Si, 

S 2 , 

S3, 

13 

98.3 

12951 

2011 Mar 28 

SI, 

S 2 , 

S3, 

12, 13 

73.7 

12952 

2011 Apr 5 

SI, 

S 2 , 

S3, 

12, 13 

142.3 

12953 

2011 Apr 7 

SI, 

S 2 , 

S3, 

12, 13 

31.7 

13246 

2011 Mar 30 

SI, 

S 2 , 

S3, 

12, 13 

45.0 

13247 

2011 Mar 31 

SI, 

S 2 , 

S3, 

12, 13 

35.7 

13253 

2011 Apr 8 

SI, 

S 2 , 

S3, 

12, 13 

116.7 

13255 

2011 Apr 10 

SI, 

S 2 , 

S3, 

12, 13 

43.0 


TABLE 2 

Properties of the SHOCKS a 


ID 

r b 

<j c 

PSF d 

\ . e 

/v in^-out 

crt f 

P2 / Pi S 

M h 

+ i 

page 


(kpc) 

(kpc) 

(kpc) 

(pc) 

(km s -1 ) 



( 10 7 yr) 

1 kpc, SE 

1.37 

0-08t^4 

0.11 

3.3 

150 ± 80 

2 06 + °' 08 
Z.OO —o 12 

l-78to.n 

0.17 

10 kpc, SE 

9.2 

0-35t“;°8 

0.15 

25 

80 ±20 

+0.04 
± - 1 °-0.03 

i e;9+0.03 
- 1 -- 0 -0.03 

1.3 

10 kpc, NW 

12.5 

0 42+ 0 - 12 

u -^ z -o.n 

0.07 

36 

70 ±20 

1 74+0.03 
‘^-0.02 

i e;9+0.03 
± - OZ -0.01 

1.8 

30 kpc, SE 

27.7 

< 1.27 

0.31 

88 

< 76 

1 26 +0 ' 08 
± -^ u —0.06 

1 1 7+0-05 
± - ±/ -0.04 

5.2 

30 kpc, NW 

29.9 

< 1.03 

0.06 

163 

< 62 

1 41+ 0 - 12 

± -^ i -0.13 

i 90 +O.O 8 

i,zo —0.09 

5.1 


a Error ranges are 90% confidence intervals. 

b Distance from the adopted center (15 b 01 m ll s .160, ±l°42 / 06 // .88) to the shock front. 

c Shock width, obtained by fitting a Gaussian smoothed density jump model to the integrated emission 
measure profile. 

d Size of the Chandra PSF in the region of the shock edge (50% EEF). 
e Mean free path of particles crossing from the post- to pre-shock region. 

f RMS speed of turbulence in the ICM as estimated from the observed shock widths and locations. 
g Density jump at shock front. 
b Mach number. 

1 Shock age, estimated as the travel time from the current position to the center point of the elliptical shock 
edge. These ages are expected to be upper limits (except for possibly the 1 kpc shock), as they assume a 
constant Mach number. 


TABLE 3 

Shock Heating and Radiative Cooling 


ID 

^COOl a 

( 10 8 yr) 

Shocks/t coo i b 

AQ/E C 

1 kpc Shock 

1.4 

7 

9 

10 kpc Shock 

9.2 

46 

21 

30 kpc Shock 

22.2 

111 

143 


a Local cooling time of the gas, just outside of the 
shock edge. 

b Number of shocks per cooling time, assuming an 
outburst repetition rate of one every 20 Myr. 
c Reciprocal of the fraction of the total thermal en¬ 
ergy E added by each shock, roughly equal to the 
number of shocks per cooling time required to offset 
radiative cooling. 
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TABLE 4 

Properties of the X-ray Cavities 


ID 

a a 

(kpc) 

b h 

(kpc) 

r c 

(kpc) 

TV d 
■C'int 

(10 56 erg) 

Inner, SW 

0.95 

0.95 

1.3 

0.5 

Inner, NE 

1.03 

0.93 

1.4 

0.6 

Middle, SW 

3.9 

3.9 

7.7 

13.1 

Middle-l e , NE 

2.9 

2.2 

4.9 

5.5 

Middle-2 e , NE 

2.8 

2.4 

9.3 

3.8 

Outer, SW 

5.2 

3.0 

22.2 

3.9 

Outer, NE 

8.0 

4.4 

18.0 

15.6 


a Semi-major axis. 
b Semi-minor axis. 
c Distance from central AGN. 

d The internal energy of the cavity, estimated as 3 PV. 
e Part of a “split” cavity. 
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Fig. 1.— Exposure corrected, background subtracted, 0.3-3 keV Chandra image of the central region of N5813, unsmoothed and with 
point sources included (1 pixel = 0.5"). 
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Fig. 2 .— Left: Exposure corrected, background subtracted, 0.3-3 keV Chandra image, with point sources removed and smoothed with 
a a = 1.5 ,/ Gaussian. The image shows bright rims surrounding an inner pair of cavities, a prominent elliptical edge surrounding a pair of 
cavities at intermediate radii (with the more obvious cavity to the SW and the NE cavity apparently broken into two connected cavities), 
and a subtle outer edge associated with a faint pair of outer cavities (with the more obvious cavity to the NE). Right: X-ray image divided 
by a 2D fitted beta model and smoothed with a cr = 6" Gaussian, shown on the same scale. The outer cavities and edges are more clearly 
seen in this residual image, while the inner cavities are not visible due to the larger smoothing scale and saturation of the color scale. The 
image also reveals a faint “channel” of decreased surface brightness extending to the north, apparently connected to the NE outer cavity. 



Fig. 3.— Left: Background and exposure corrected 0.4-7.2 keV XMM-Newton image of N5813. Right: Smoothed Chandra image shown 
on the same scale, with the intensity scale chosen to better show the faint, outer emission. 
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Fig. 4.— XMM-Newton image shown in Figure [3] smoothed and with the intensity scale chosen to show the faint emission at large 
radii beyond the Chandra FOV. There are no additional cavities or shock edges visible beyond the inner features identified in the Chandra 
images. 
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Fig. 5.— Smoothed (left) and contour binned (right) temperature maps. The locations of the 10 kpc and 30 kpc shock fronts are 
indicated with dashed black lines, and the cavity locations with black ellipses. Both maps clearly show temperature increases associated 
with the 1 kpc and 10 kpc shocks, and hint at increases associated with the 30 kpc shock (particularly to the NW). Also visible is a SW 
to NE plume of cool gas that has been uplifted by the buoyantly rising cavities. 



Fig. 6. — Left: High-resolution smoothed temperature map of the core region. Right: The contour binning temperature map shown on 
the same scale. The 10 kpc shock fronts are indicated as in Figure [5] 
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Fig. 7.— Pseudo-pressure (left) and pseudo-entropy (right) maps, in arbitrary units. The pressure map was calculated as kTA 1 / 2 and 
the entropy map as kTA ~ 1 / 3 , where A is the APEC normalization scaled by the area of the extraction region. Shocks and cavities are 
indicated as in Figure [5] The pressure jumps are visible at the 1 kpc and 10 kpc shock fronts. There are no visible entropy jumps at 
the shock fronts, consistent with the expectation that entropy jump amplitudes are small for weak shocks compared with those of the 
temperature and pressure jumps. 



Fig. 8 .— Pseudo-pressure (left) and pseudo-entropy (right) maps corresponding to the smoothed temperature map of the core shown in 
Figure [m (with the same regions overlaid), created as in Figure |7| Pressure increases are clearly seen at the 10 kpc shock fronts and in the 
bright, shock heated rims surrounding the inner cavity pair. There is some structure visible in the core of the pseudo-entropy map (seen 
in blue as filament-like structures), likely resulting from central gas that has been pushed out and uplifted from the core by the expanding 
and rising X-ray cavities. 
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Fig. 9.— Core temperature (left) and pseudo-entropy (right) maps with Ha contours from Rll overlaid. Units are keV for temperature 
and arbitrary for entropy. The SE 10 kpc shock is indicated by the dashed line. The Ha emission follows the cool plume of gas up to 
the intermediate cavities, and traces the SW inner cavity. There is a correlation with the detailed, filamentary structure of the central 
low-entropy gas in the pseudo-entropy map (e.g., with the low-entropy filament just NE of the core). 
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Fig. 10.— Smoothed abundance map, with cavity and shock regions overlaid as in Figure [5] The apparent low abundance in the plume 
of uplifted cool gas (extending NE and SW of the core) is an artifact of fitting a single temperature model to multi-temperature spectra 
(see text). 
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Fig. 11.— Azimuthally averaged radial profiles for (top to bottom) temperature, electron density, pressure, and entropy, extracted from 
circular annuli. The top panel shows the projected (black circles) and deprojected (red triangles) temperature profiles. All other panels 
show deprojected values. The vertical lines mark the positions of the shock fronts to the SE (dotted) and NW (dashed) determined by 
fits to the surface brightness profiles in sectors (see §(5f. Increases in temperature associated with each pair of shock fronts (i.e., each full 
elliptical edge) are clearly seen even in the azimuthally averaged, projected temperature profile. The “kink” in the deprojected profiles is 
likely due to the breakdown of the assumption of spherical symmetry at the location of the bright, sharp shock edges at ~ 10 kpc. 
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Projected temperature profiles across the NW (left) and SE (right) 10 kpc 
positions. 
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Fig. 14.— Integrated emission measure profile across the NW middle 10 kpc) shock. The black solid line shows the fitted projected 
model consisting of a discontinuous power-law density profile (with separate inner and outer power-law slopes). The red dashed line shows 
the fit with an identical model smoothed with a Gaussian, with the Gaussian width allowed to vary. The vertical dotted line shows the 
location of the shock front. The Gaussian smoothed model gives a better fit to the data. 
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Fig. 15.— 0.3-3.0 keV surface brightness profile in azimuthal bins across the northern channel of decreased surface brightness shown in 
Figure |2] The angle © is measured from east to north. There is a clear dip at the location of the northern channel, indicating that this is 
a real feature. 






















